void SAMPLE(I, En, Vir)
{
	int I;
	double En, enp, Vir, press, CORP, vol, rho;
	if (NPART != 0)
	{
		enp = En/double(NPART);
		vol = BOX * BOX * BOX;
		pres = (NPART/vol)/BETA + Vir/(3*vol);
		rho = NPART/vol;
		if (TAILCO)
			press = press + CORP(RC, rho);
	}
	else{
		enp = 0;
		press = 0;
	}
	printf(I, enp, press);
}
